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The explosion energy of early stellar populations: The Fe-peak 
element ratios in low metallicity damped Lya systems* 
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ABSTRACT 

The relative abundances of the Fe-peak elements (Ti-Zn) at the lowest metallicities are in- 
timately linked to the physics of core-collapse supernova explosions. With a sample of 25 
very metal-poor damped Lya systems, we investigate the trends of the Fe-peak element ratios 
with metallicity. For nine of the 25 DLAs, a direct measurement (or useful upper limit) of 
one or more of the Ti,Cr,Co,Ni,Zn/Fe abundance ratios could be determined from detected 
absorption lines. For the remaining systems (without detections), we devised a new form of 
spectral stacking to estimate the typical Fe-peak element ratios of the DLA population in this 
metallicity regime. We compare these data to analogous measurements in metal-poor stars 
of the Galactic halo and to detailed calculations of explosive nucleosynthesis in metal-free 
stars. We conclude that most of the DLAs in our sample were enriched by stars that released 
an energy of < 3.0 x 10 51 erg when they exploded as core-collapse supernovae. We find one 
DLA that appears to be modestly supersolar in Ni/Fe, and tentatively suggest that this system 
may have been enriched by hypernovae. We conclude by discussing the exciting prospect of 
measuring Fe-peak element ratios in damped Lya systems with Fe/H < 1/1000 of solar when 
30-m class telescopes become available. Only then will we be able to pin down the energy 
that was released by the supernovae of the first stars. 

Key words: galaxies: abundances - galaxies: evolution - quasars: absorption lines - super- 
novae: general 



1 INTRODUCTION 

The physical mechanism that drives a core-collapse supernova 
(CCSN) explosion is a challenging mystery that remains un- 
solved. However, it has been widely recognised for some time 
that the relative abundances of the Fe-peak elements (Ti-Zn), 
which predominantly originate from explosive nucleosynthesis 
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(see e.g. IWooslev & WeaverJll995h . can provide unique insights 
into the inner workings of a massive star in its death throes 
dUmeda & Nomotol l2002h . In particular, the ratios of the Fe- 
peak element yields depend sensitively on the progenitor mass, 
the explosion energy, and the amount of material that falls back 
onto the remnant (i.e. the mass-cut). The uncertain physics of 
the explosion mechanism forces model nucleosynthesis calcu- 
lations to parameterise the explosion. There are two favoured 
prescriptions to define the central engine of the explosion: the 
first mecha nism involves depos i ting momentum at a given mass 
coordinate dWooslev & Weaverj 1 19951 ; IChieffi & Limongil |2004 
iHeger & Wooslevl l20ld)T whereas the s econd is a deposition of 
thermal energy dNakamura et al. I l200ll) . Both prescriptions ap- 
pear to give consistent yields for the lighter a-elements, although 
they produce quite different abundances o f the Fe-peak elements 
dAufde rheide. Ba ron. & Thielemannlll99ll) . This is not surprising, 
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given that the Fe-peak elements are intimately linked to the details 
of the explosion. 

It is reassuring that one of these prescriptions (momentum 
deposition) is able to do a remarkably good job at reproducing 
the large variety of abundance patterns that massive Population 
III stars (10-100M G ) encoded in the second generation of stars 
(see, for example, the comparisons performed bv lHeger & Wooslevl 
l20ld) . These second gener ation stars have largely been discov ered 
through the dedicated HK (Beers, Preston, & Shectman 1992) and 
Hamburg/ESO (HES; IChristlieb et alj|200ll) surv eys, with detailed 
abund a nce analyses being n otabl y conducted by | Mc William et al .1 
dl995ri , ICavrel et al.l d2004l) . and lLai et alj J2008h and references 
therein. The observations have uncovered some unexpected trends 
in the relative abundance of Fe-peak elements with decreasing 
metallicity, which we now briefly discuss. 

The first detailed study of the Fe-peak elements in stars 
with met allicities < 1/100 of solar ([Fe/H] < -2) was con- 
ducted bv lMcWilliam et all d 1 9951) . Their sample comprised the 33 
most metal-poor halo stars known at the time, extending down to 

[Fe/H] 4.0. These authors noted changes in the general trends of 

some of the Fe-peak element ratios for metallicities [Fe/H] < -2.4. 
In particular, they found that the Cr/Fe ratio, which is roughly so- 
lar in stars with [Fe/H] > -2.4, decreases to [Cr/Fe] ^ -0.5 at 
[Fe/H] = -4.0, albeit with so me scatter. This tren d was later con- 
firmed bv lCavrel et all j2004) and lLai et alj J2008h . whose data ex- 
hibited reduced scatter. Unfortunately, there is still a worrisome 
disagreement between the Cr/Fe ratios deduced from either the 
Cri or Cm lines; differences of up to 0.3 dex for metallicitie s 
[Fe/H] < -2.4 may be expected dSobeck. Lawler. & Sneder]|2007h . 
and the true Cr/Fe relative abundance in the low metallicity regime 
may in fact not be significantly different from the solar value. 

The Ni/Fe ratio, on the other hand, exhibits no discernible 
trend away from the solar value down to the lowest metallici- 
ties probed by the above surveys. This is quite unlike the Zn/Fe 
ratio, which is roughly solar in the metallicity interval 0.0 > 
[Fe/H] > -2.0, but then gradually increases wit h decreasing metal- 
licity reaching [Zn /Fe] a +0.7 at [Fe/H] = -4.0 dPrimas et al.l200d ; 
ICavrel et al.ll2004l) . If such a behaviour were due to CCSN nucle- 
osynthesis, the excess p roduction of zinc is likely the result of a 
strong a-rich freeze out dWooslev & WeaveJI 19951), which in turn 
is linked to the explosion energy. In fact, lUmeda & Nomotd llwA) 
invoked the dominance of hypernovae — a CCSN explosion with an 
order-of-magnitude more energy released than normal — at the low- 
est metallicities to explain the elevated [Zn/Fe] values. This inter- 
pretation is also consistent with the observed rise in the [Co/Fe] val- 
ues exhibited by the same group of stars. Despite these significant 
advances in the stellar abundance measurements, it would be highly 
advantageous to have complementary probes to independently con- 
firm the trends uncovered in the most metal-poor stars. This will 
ensure that we are not biased by systematics that may affect the 
modelling of the absorption lines in the stellar atmospheres. 

With this goal in mind, we and other groups have re- 
cently focussed our attention on the most metal-poor clouds 
of neutral ga s at redshifts z ~ 2 -3 jPettini et alj 1 20081; 



Penprase et al] 120081 ; ICooke et alj l2011al lbl; see also iBecker et alj 
20121) . Such clouds are identified amongst the damped Lya sys- 



first ions (e.g. IVladilo et alj 1200 ll) . This simple ionization struc- 
ture makes the measurement of element abundances straightfor- 
ward. I n a recently comple ted survey of 22 DLAs with [Fe/H] 
< -2.0. ICooke et alj d201 lbl) noted a good agreement in the abun- 
dances of C, O and Fe between halo stars and DLAs of comparable 
metallicity. This is perhaps not surprising, since high redshift very 
metal-poor (VMP) DLAs likely contain the gas from which the 
most metal-poor stars m ay have later condensed. As discussed by 
IPettini & Cookd d2013l) , VMP DLAs are therefore an ideal, com- 
plementary probe to measure the relative abundance of the Fe-peak 
elements, and compare them to those found in Galactic halo stars. 

This is the subject of the present paper, where we attempt to 
measure Fe-peak element ratios in a sli ghtly expanded vers ion of 
the sample of VMP DLAs considered bv lCooke et"al]d2011bl) . This 
is a difficult prospect because most Fe-peak elements are less abun- 
dant than Fe by one to two orders of magnitude and the absorption 
lines of their dominant ion stages in VMP DLAs are generally too 
weak to be detected with current observational facilities. To cir- 
cumvent this issue, we introduce a form of spectral stacking that 
allows us to measure typical Fe-peak element ratios for the whole 
sample of VMP DLAs. In Section 2, we briefly describe the sam- 
ple. In Section[3]we introduce the spectral stacking technique and 
present the results of our analysis. We then compare in Section [4] 
the DLA Fe-peak element ratios with those measured in metal-poor 
halo stars. Finally, in Section [5] we compare these ratios with the- 
oretical yields calculated with the current generation of explosive 
nucleosynthesis models for metal-free stars, before concluding in 
Section|6] 



2 THE SAMPLE 

The majority of the VMP DL As considered in the present work 
were drawn from the survey bv lCookeetal] d2011bl) which com- 
prises 22 DLAs with metallicities less than 1/100 of solar. We ex- 
cluded three DLAs from that sample. Two of these were excluded 
because no Fen absorption lines were detected. In the third case, 
we found that the wavelength calibration of the spectrum was too 
uncertain to reliably stack portions of the spectrum for the pur- 
pose of detecting a signal from weak absorption lin es (see Sec- 
tion 13. U . We augmented the VMP DLA sample of ICooke et alj 
( 12011b) with two further systems from our own work, one pub- 
lished ( C ooke. Pettini. & Murphy|[2012T) and one which is the sub- 
ject of a forthcoming study (Cooke et al., in preparation). Further- 
more, we have also included five VMP D LAs reported in the li tera- 
ture that didn't meet the selection cuts of lCooke et alj d2011bl) . but 
show clear detections of absorption lines from Fe-peak elements. 
Thus the total sample considered here consists of 25 VMP DLAs. 
Details of these systems are collected in Table Q] 

The composition of the DLAs is derived by searching for 
the absorption lines of elements in their dominant ionization state 
(usually the neutral or singly ionized species for DLAs) against 
the emission spectrum of a more distant, background QSO. In all 
cases, the QSOs were observed with either the European Southern 
Observ atory's (ESQ) Ultra violet and Visual Echelle Spectrograph 



terns (DLAs), absorbers with neutral h ydrog en column densities 
N(Hi) > 10 203 atoms cm' 2 (IWolfe et al|l98d ; for a review of this 
topic see I Wolfe. Gawiser. & Prochaskall2005l) . At such large col- 
umn densities the clouds are self-shielded from the metagalac- 
tic ionizing background, and most elements are concentrated in 
a single dominant ionization state, usually the neutrals or the 



(UVES. lDekker et al.l20od) . or the W. M. Keck Observatory 's High 
Resolution Echelle Spectrometer (HIRES. IVogt et all 19941) . These 
spectra typically have resolving power R ~ 40 000, and signal-to- 
noise ratio S/N ~ 30 per pixel (~ 2.5 k m s" 1 ). Detail s of the data 
acquisition and reduction can be found in ICookeetalJd2011bh . 

For the present study, we extracted a +300 km s"" 1 window 
around each of the absorption lines of interest and fitted a low or- 
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Table 1. The DLA sample 



3.1 Stacking Technique 



QSO Name 


- m 


Zabs 


N(Fe n) 
8 cm- 2 


, N(Hl) 

° cm z 


S/N a 


J0035-0918 


2.420 


2.3400972 


12.99 


20.55 


17 


Q01 12-306 


2.985 


2.4184400 


13.33 


20.50 


23 


JO 140-0839 


3.716 


3.6966000 


12.77 


20.75 


31 


1083 1+3358 


2.427 


2.3036420 


13.33 


20.25 


15 


00913+072 


2.785 


2.6183762 


12.99 


20.34 


47 


J 100 1+0343 


3.198 


3.0784130 


12.50 


20.21 


40 


J1037+0139 


3.059 


2.7048700 


13.53 


20.50 


45 


Ql 108-077 


3.922 


3.6076700 


13.88 


20.37 


49 


Jl 1 12+1333 


2.420 


2.2709396 


13.58 


20.44 


53 


J1337+3152 


3.174 


3.1676670 


13.14 


20.41 


19 


J 1340+ 1106 


2.914 


2 5079179 


13.49 


20.09 


50 


J1340+1106 


2.914 


2.7958257 


14.32 


21.00 


50 


J1358+6522 


3.173 


3.0672953 


13.10 


20.47 


30 


J1419+0829 


3.030 


3.0498400 


13.54 


20.391 


43 


J1419+5923 


2.321 


2.2472964 


14.35 


20.95 


18 


J1558-0031 


2.823 


2.7023887 


14.11 


20.67 


31 


J1558+4053 


2.635 


2.5537663 


13.07 


20.30 


IS 


J 1944+7705 


2.994 


2.8443000 


13.24 


20.27 


110 


Q2059-360 


3.090 


3.0828330 


14.48 


20.98 


27 


J2155+1358 


4.256 


4.2122097 


12.93 


19.61 


22 


Q2206-199 


2.559 


2.0762272 


13.33 


20.43 


100 



''Indicative signal-to-noise ratio at 5000 A (or 6000 A in the case of 
J0140-0839, Ql 108-077, J1337+3152, J1340+1106 and J2155+1358). 



der polynomial to regions deemed to be free of absorption. The 
details of this continuum fit (i.e. the polynomial order and pixels 
used in the continuum fit) are stored, and are used later to derive 
the systematic errors introduced by our choice of continuum. The 
absorption lines targeted are listed in Table0 vacu um wavelengths 
and oscillator st rengths are from t he com pilation bv lMortoril d2003h 
with updates by Jenk ins & Tripd d2006l) . 



3 ANALYSIS AND RESULTS 

The primary goal of this study is to measure the relative abundances 
of a few Fe-peak elements with atomic number in the range 22 (Ti) 
to 30 (Zn), in the lowest metallicity systems known at redshifts z <; 
2. Of these elements, we adopt Fe as the 'reference' element since 
it is the most abundant of the Fe-peak elements, and it is detected 
in all of the DLAs in our sample with at least 5cr confidence. 

The remaining Fe-peak elements are all less abundant than 
Fe, and their absorption lines are rarely seen in the most metal- 
poor DLAs; the only reported detections in individual VMP DLAs 
are those discussed here. The reason for such dearth of data can be 
appreciated by considering the following. Even the strongest Ni n 
absorption line in Table|2] Niii/11317.217, would require a signal- 
to-noise ratio S/N > 110 per resolution element to be detected at 
the 5cr level in a DLA with [Fe/H] < -2 and log N(Hi) = 20.3. 
Such a high S/N ratio is beyond current observational capabilities 
given that most VMP DLAs are found in the spectra of faint QSOs 
( V > 18). Consequently, all measures of [Ni/Fe] in individual VMP 
DLAs reported to date refer to higher (and rarer) column density 
systems than the DLA threshold log JV(Hi) = 20.3. 

To circumvent this issue, we developed a new type of stacking 
technique to measure, or place limits on, the ratios of several Fe- 
peak elements in the DLAs in our sample. 



The stacking technique that we now describe can be used to com- 
bine either: (i) multiple transitions of a given ion for a single DLA; 
or (ii) all transitions of a given ion in all DLAs (assuming that 
Fe-peak element ratios do not vary significantly between different 
VMP DLAs). In individual DLAs where we detect at least one ab- 
sorption line from an Fe-peak element (other than Fe), we used the 
former approach and stacked all transitions of that ion covered by 
our spectrum of the DLA, whether detected or not. This has the 
effect of utilising all of the information available to improve the 
precision of the ion column density determination over that achiev- 
able from the analysis of only the absorption lines that have been 
detected. For DLAs where a given ion has not been detected, we 
adopted the latter approach and produced a 'population' stack. 

One advantage of working with weak, undetected absorption 
lines is that they fall on the linear part of the curve of growth^. In 
this regime, one can simply relate the rest-frame equivalent width 
of an absorption line to the column density (cnr 2 ) of element X in 
ionization state n by: 



EW(Xn) 



N(XN)fA 2 



(1) 



1.13 x 10 20 ' 

where / and A are the oscillator strength and the rest wavelength 
of the atomic transition (in A) respectively; values relevant to this 
work are collected in Table [2] This technique has the clear advan- 
tage of being independent of the cloud modelling, and the spectral 
resolution. Given that in all cases we have least one unsaturated 
Fen absorption line from which the column density of singly ion- 
ized iron, A'(Fen), can be deduced with confidence, equation {TJ 
can be rearranged as follows: 



( N(Xu) \ 
U(Fen)l 



EW pp (Xn) 



1.13 x 10 20 
A'(Fen)//! 2 ' 



(2) 



where EW pp is the equivalent width per pixel. Thus, for each pixel 
in the spectrum where we see absorption by Xn, the rest-frame 
equivalent width of this pixel alone (weighted by the fraction on 
the right side of equation represents the fractional contribution 
to the total column density ratio of species X n relative to Fe n. If the 
width of each pixel is measured in km s -1 , the quantity on the left 
side of equation has the units knT 1 s. By integrating equation 
over all pixels where the absorption takes place then gives the final 
value of N(K N)/AT(Fe n). 

For a given DLA, equation can be used to create a 'stacked' 
spectrum for all lines of a given species. As an example, the suite 
of Ni ii lines available in many of the DLA spectra can be combined 
to provide a single, more precise measure of the Ni n/Fe n ratio. To 
do this, we first convert the spectrum in the vicinity of each tran- 
sition to a common velocity scale. The zero point of the velocity 
scale is established from the profiles of the Fe n lines available. We 
then calculate the rest-frame equivalent width for each pixel for all 
transitions and use equation to convert this to a column density 
ratio. Finally, we determine the weighted mean of all transitions for 
a given species, where the weights correspond to the inverse vari- 
ance of each pixel. In this way, we ensure that the final stack has 
the highest achievable S/N. We stress that this technique can be im- 
plemented for a single DLA without any assumptions. Therefore, 
for DLAs where we detect at least one transition of a given Fe-peak 
element (apart from Fe itself), we will use this technique to stack 

1 This statement does not necessarily apply to low resolution spectra, but 
is valid at the high spectral resolution of our echelle observations. 
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Table 2. Adopted metal line laboratory wavelengths and oscillator strengths 11 



Ion 


Wavelength 


/ 


Ion 


Wavelength 


/ 


Ion 


Wavelength 


f 




(A) 






(A) 






(A) 




Si ii 


1808.01288 


0.00208 


Con 


1424.7866 


0.0109 


Nin 


1370.132 


0.05 


Sn 


1250.578 


0.00543 


Con 


1466.2110 


0.031 


Nin 


1393.324 


0.0101 


Sn 


1253.805 


0.0109 


Con 


1480.9546 


0.0119 


Nin 


1454.842 


0.0323 


Sn 


1259.5180 


0.0166 


Con 


1552.7624 


0.0116 


Nin 


1502.1480 


0.0133 


Ti ii 


1910.9538 


0.098 


Con 


1572.6480 


0.0120 


Nin 


1709.6042 


0.0324 


Tin 


1910.6123 


0.104 


Con 


1574.5508 


0.025 


Nin 


1741.5531 


0.0427 


Cm 


2056.25693 


0.1030 


Con 


1941.2852 


0.034 


Nin 


1751.9157 


0.027 


Cm 


2062.23610 


0.0759 


Con 


2012.1664 


0.0368 


Znn 


2026.13709 


0.501 


Cm 


2066.16403 


0.0512 


Nin 


1317.217 


0.057 


Znn 


2062.66045 


0.246 



"From Morton (2003) with updates by Jenkins & Tripp (2006). 



all the available transitions in that DLA spectrum and obtain the 
value of A'(XN)/A'(Feii) appropriate to that DLA. 

The stacking method just described can only be applied to 
the few cases where we do detect at least one absorption line of 
a given Fe-peak element. However, for most VMP DLAs in our 
survey, even stacking all the (undetected) absorption lines of a 
given ion in the spectrum will not improve the S/N sufficiently to 
lead to a meaningful measure of yV(XN)/W(Fen). In these circum- 
stances, however, we can still apply the above method to construct 
a stacked spectrum of all lines of a given ion for all DLAs with 
non-detections. In this way, we can derive an averaged measure of 
A'(XN)/A'(Feii) (or an upper limit to this ratio) for the population 
of VMP DLAs in our sample. Such an average is still instructive, 
given the relatively small scatter in the ratios of most Fe-peak el- 
ements exhibited by Galactic halo stars of comparab le metallicity 
dCavrel et alj|2004l ; lLai et a"i]|2008l: lYong et al.ll2012» . 

There are two potential systematics that may introduce a bias 
when stacking the absorption lines for a selection of DLAs with a 
variety of metallicities and Hi column densities. The first concern 
is that certain species may be depleted (relative to Fe) onto dust 
grains, an effect which is known to depend on metallicity. Fortu- 
nately, when the DL A metallicity is [Fe/H] < -2.0, dust deple tions 
become unimportant jPettini et alJl997l ; lAkerman et al.l2005t) . The 
other concern is that the radiation field to which the gas is ex- 
posed may vary from DLA to DLA, depending among other factors 
on the Hi column density. However, for the ions of the elements 
under consideration here, which are the dominant ions in neutral 
gas, the differential ionization corrections required to convert the 
measured iV(XN)/W(Fen) to the element abun dance ratio X/Fe are 
minimal and can be safely neglected (see e.g. IVladilo et alj|200ll ; 
ICooke et alfeoioh . 

Thus, in what follows we assume N(X N)//V(Fe ii) = 
N(X)/N(Fe). Comparing the above ratio to the solar value then 
gives the quantity of in terest, [X/Fe] = log (jV(X)//V(Fe)) - 
log(X/Fe) Q ). We adopt the lAsplund et al.l J2009I) solar abundance 
scale, taking either the photosphe ric, meteoritic, or an average o f 
the two determinations following lLodders. Plame. & Gaill d2009l) . 
The solar abundances adopted in this work are listed in Table[3] 

The most likely source of systematic error is in the choice of 
the continuum around each absorption line, which is especially im- 
portant when dealing with weak lines. To estimate the additional 
uncertainty (and bias) that is introduced by our placement of the 
continuum le vel, we adopted a Monte Carlo approach similar to that 
described by IPettini & Cookej d2012h . We used the measured line 
equivalent widths (or upper limits) to generate a series of line pro- 



Table 3. Adopted solar abundances 11 



X 


log(X/H) 


X 


log(X/H) 


Si 


-4.49 


Fe 


-4.53 


S 


-4.86 


Co 


-7.07 


Ti 


-7.09 


Ni 


-5.79 


Cr 


-6.36 


Zn 


-7.37 



a From Asplund et al. (2009). 

files for each of the transitions used in the stack. We then combined 
these profiles with the stacking technique described above — first 
without any continuum adjustments and then again with the identi- 
cal continuum adjustments applied to the real data (see Sec.0. This 
exercise was repeated 1000 times, and the deviations between the 
equivalent widths derived with and without the continuum adjust- 
ments were taken to be a measure of the systematic uncertainty. In 
general, we found that the systematic errors are of the same order as 
the random errors. Perhaps more importantly, we also verified that 
in all of the cases described below there were no significant system- 
atic biases introduced by our estimate of the continuum level. 



3.2 J1340+ 1 106 - A test case 

We tested the validity of this technique by applying it to one of the 
few clear detections of Cr n and Ni n absorption in our sample. With 
a relatively high column density A'(Hi) = 1 x 10 21 cm" 2 , the DLA 
at redshift z a bs = 2.79583 toward J1340+1106 shows an unusually 
high number of met al lines for its low metallicity [Fe/H] = -2.15 
dCookeetal.l2011bl) . From the ana lysis of 36 absorption lines from 
1 1 ions dominant in H i regions, ICooke et alj J201 lbh derived a 
'cloud model' which successfully reproduced the line profiles. This 
model was used to derive the ion column densities, including those 
of Cm and Ni n. 

In Fig.Q] we have reproduced the stacks of Cr n and Ni n lines 
derived from equation 10, showing clear detections near relative 
velocity v = km s _1 . The total Crii/Feii and Nin/Fen ratios are 
obtained by integrating these profiles over the velocity range of the 
absorption, defined by the extent of the stronger absorption lines 
of Fen. We find [Cr/Fe] = 0.08 ±0.11 and [Ni/Fe] = 0.03 ± 0.03. 
These values compare well with [Cr/Fe] = 0. 13±0.11 and [Ni/Fe] = 
0.02 ± 0.03 deduced bv lCooke et all d2011bl) from profile fitting of 
the absorption lines. 
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J1340+1106, z abs =2. 79583 



£ 

CO 

o 

X 

\ 
u 




■100 -50 50 100 
Velocity (km s -1 ) 



J1340+1106, z abs =2. 79583 




e 

CO 

o 

X 
CD 



-100 -50 50 100 
Velocity (km s -1 ) 

Figure 1. The stacked spectrum (black histogram) that is used to measure 
the Cr/Fe and Ni/Fe ratios (top and bottom panel respectively) in the DLA 
at ; abs = 2.79583 towards J1340+1106. The shaded red bands show the 
±1lt uncertainty, whereas the long-dashed blue lines indicate the levels for 
no Cr ii and Ni n absorption. The black 'X' at +55 km s~ 1 in the top panel 
marks an unrelated feature. 



3.3 Individual systems with detections 

We now apply the stacking technique described in Section [37X1 to 
the five DLAs in our sample that: (i) either show absorption by at 
least one Fe-peak element (other than Fe), or (ii) where the data are 
of sufficient quality to allow a useful upper limit to be set. For this 
to be the case, these VMP DLAs either have a large H i column den- 
sity, or their spectra have been recorded at an unusually high S/N, 
or both. Even in these most favourable cases, however, the only Fe- 
peak elements detected are Ni and Cr, the two most abundant after 
Fe. In Table|4]we list five detections (four for Ni n and one for Cr n) 
and three upper limits (one for Nin and two for Cm). For each 
system, we give [Fe/H], the significance of detections of Nin and 
Cr n, and the number of transitions used in the stacked spectrum. In 
all cases, we integrated the stacked profile over the velocity range 
exhibited by the stronger Fe n absorption lines, which is provided 
in column 6 of Table|4] 

The stacked profiles for these five DLAs are shown in Fig. [TJ 
and Fig. [2] All of the right panels in Fig. [2] have significant fea- 
tures (see penultimate column of Table 0). In the lower-left and 
middle-left panels of the figure, we see no discernible absorption 
and deduce corresponding 3<x upper limits to the ratios Ni/Fe and 



Cr/Fe. The top-left panel shows a feature that is 'detected' at the 
3cr level from which we conservatively derive a 3<x upper limit to 
[Cr/Fe]. 



3.4 Systems with non-detections 

The remaining 16 DLAs in our sample do not show individual ab- 
sorption lines from any Fe-peak element besides Fe; these systems 
are now used to produce a mean 'population' stack. As discussed 
earlier, in generating this stack we assume that VMP DLAs do not 
exhibit a wide range of [X/Fe] values over the metallicity range 
probed here, -3.5 < [Fe/H] < -2.0. Cases with individual detec- 
tions described in Section [331 were not included in the stacks. By 
inspecting the Fe n absorption lines for all systems in the stack, we 
determined that most absorption occurs over the velocity interval 
Av = +25 km s from the reference redshift. We then confirmed 
that all of the undetected absorption lines used in the stack con- 
tained no significant absorption within this velocity interval (i.e. 
all pixels are within +2cr of the continuum level in this velocity 
interval). Furthermore, to be certain that we are not including spu- 
rious contributions from coincident, unrelated absorption, we re- 
jected the high and low values for each pixel in the stack. This 
procedure ensures that we suffer minimal contamination. 

The results of these stacks are reproduced in Fig. [3] Due to the 
wider variety of line profile shapes that have been used in the stack 
(which for some cases contain multiple components and extend to 
+25 km s -1 ), the features appear more broadened than those shown 
in Fig.|2]which typically exhibit a single dominant absorption com- 
ponent. The corresponding values of [X/Fe] are listed in the fourth 
column of Table|4] The metallicity [Fe/H] to which these values re- 
fer is not the same for every element X considered, because not all 
of the 16 DLAs contribute the same number of absorption lines to 
each stack. The values of [Fe/H] and their ranges listed in the fifth 
column of Table |4]under the heading 'Stacked systems with non- 
detections' were determined using a weighted bootstrap method on 
the array of Fe/H values for each stack, with /V(Fen) x(S/N) 2 as 
the weight on each line used in the stack. 

In order to check the consistency of our results, we also 
produced two stacks of lines from o--capture elements, one for 
the weak Sin A 1808 line, and the other for the Sn triplet at 
AA1250, 1253, 1259. As can be seen from Table H we find 
[Si/Fe] = +0.31 ± 0.04 ± 0.05 (random and systematic error respec- 
tively) and [S/Fe] = +0.25 ± 0.05 + 0.06. The modest enhancement 
by a factor of ~ 2 of these two cr-capture elements relative to Fe is 
co nsistent with the mea n value [(Si/Fe)] = +0.32 ± 0.09 deduced 
by |Cookeet"aTJ ( l2011bl) in the same metallicity interval. 



3.5 Other systems in the literature 

To our knowledge, this is the first detailed investigation into the 
Fe-peak element ratios in very metal-poor DLAs. This is perhaps 
not surprising, given that there are just a few VMP DLAs that 
we could find in the literature with column density measurements 
of Fe-peak elements: (1) the DLA w ith /V(Hi)= 21.41 + 0.08 at 
z abs = 3.3901 towards Q0000-2619 dMolaro et al1l2000l) : (2) the 
z abs = 3.915 DLA tow ards Q0255+00 with /V(Hi)= 21.30 ± 0.05 
jProchaska et alj200ll) : (3) the sys tem at z abs = 2.615 with N(R i)= 
21.30 ± 0.10 towards Q2348-01 jProchaska etai]|200l[) : and (4) 
the DLA to wards Q2348-1444 at z abs = 2 .279 with N(Hi)= 
20.59 + 0.08 jPessauges-Zavadskv et~aT1l2007t) . Other relevant de- 
tails for these absorbers are collected in Table|4] under the heading 
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Figure 2. Stacked spectra (black histograms) for individual DLAs in which at least one absorption line from an Fe-peak element (other than Fe) was detected, 
or a useful upper limit could be determined. The red shaded regions correspond to the ±l<x random errors. The three panels on the right show significant 
detections of Ni n, while from the three panels on the left 3<x upper limits to [Cr/Fe] and [Ni/Fe] can be determined. Black 'X' symbols denote unrelated 
absorption features. 



'Literature sample'. All literature values have been adjusted to the 
same solar abundance scale adopted here. 



4 COMPARISON WITH LOW-METALLICITY STARS 

Our knowledge of the Fe-peak elements in the VMP regime 
([Fe/H] < -2.0) is currently based primarily on observations of the 
most metal-poor stars in the halo of our Galaxy and its companions. 
VMP DLAs offer a new probe to verify independently the validity 
of the stellar trends with metallicity. One advantage of using DLAs 
to measure relative element abundances is that the derived ratios 



are model independent; in metal-poor stars, systematic uncertain- 
ties in the abundances may arise for some absorption lines that are 
modelled in one-dimension, or if the assumption of local thermo- 
dynamic equilibrium (LTE) is not valid. 

To obtain independent confirmations of the stellar trends, 
we compare in Fig. [4] the Fe-peak ratios Ti, Cr, Ni, Zn/Fe mea- 
sured in our sample o f VMP DLA s with the values found i n 
the halo star s u rveys bv lGratton et~ai1 j2003l) . ICavrel et all ( 120040 . 
and lLai et al. I J2008I) , plott ed as open ci r cles, s quares, and di- 
amonds respectively. Since iGratton et all d2003l) did not report 
[Co/Fe] measurements, we use the recen t literature compilation by 
iBergemann. Pickering. & Gehrenl jiolch . and display these values 
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Figure 3. Stacked spectra (black histograms) for DLAs without individual detections of Fe-peak elements. The red shaded regions correspond to the ±lcr 
random errors. We have clear detections for Si/Fe, S/Fe, Cr/Fe, and Ni/Fe and 3<r upper limits for Ti/Fe, Co/Fe, and Zn/Fe. Black 'X' symbols denote unrelated 
absorption features. 
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Table 4. Summary of stacked abundances 



Name z abs X [X/Fe] a [Fe/H] Av Confidence Number 

kms~' of lines 



Individual systems with detections 



J1112 + 1333 


2.27094 


Cr 


< +0.12 


-2.33 ± 0.06 


±10 


3<r 


2 


J1112+ 1333 


2.27094 


Ni 


-0.19 ±0.12 ±0.06 


-2.33 ± 0.06 


±10 


4<r 


6 


J1340+ 1106 


2.50792 


Ni 


s; +0.02 


-2.07 ± 0.05 


+ 15 

-30 


3<r 


5 


J1340+ 1106 


2.79583 


Cr 


+0.08 ±0.11 ±0.10 


-2.15 ±0.06 


+ 13 


4 ( r 


2 


J1340+ 1106 


2.79583 


Ni 


+0.03 ± 0.03 ± 0.02 


-2.15 ±0.06 


;H 

-30 


17cr 


7 


J1419 + 5923 


2.24730 


Ni 


-0.04 ± 0.05 ± 0.03 


-2.07 ±0.10 


±20 


lOcr 


7 


J1944 + 7705 


2.84430 


Cr 


< +0.16 


-2.50 ± 0.06 


±15 


3<r 


2 


J 1944 + 7705 


2.84430 


Ni 


+0.25 ± 0.05 ± 0.03 


-2.50 ± 0.06 


±15 


9<r 


7 



Stacked systems with non-detections 

STACK Si +0.31 ±0.04 ±0.05 -2.33+!! !' ±25 13cr 6 

STACK S +0.25 ± 0.05 ± 0.06 -2.55+™ ±25 6cr 14 

STACK Ti < +0.16 -2.19!;™ ±25 3cr 16 

STACK Cr +0.1 1 ± 0.08 ± 0.09 -2.27^ ±25 5cr 8 

STACK Co < +0.53 -2.17+™ ±25 3cr 85 

STACK Ni -0.03 ± 0.04 ± 0.03 -2.19™° ±25 12cr 52 

STACK Zn < -0.04 -2.32+fjjj ±25 3cr 11 



Literature sample 



Q0000 


-2619 b 


3 


39010 


Cr 


+0.05 ± 0.04 


Q0000 


- 2619 b 


3 


39010 


Ni 


-0.22 ± 0.04 


Q0000 


- 2619 b 


3 


39010 


Zn 


-0.02 ± 0.06 


Q0255 


+ 00 c 


3 


91500 


Ni 


-0.22 ±0.10 


Q2348 


-01 c 


2 


61500 


Cr 


-0.07 ±0.11 


Q2348 


-01 c 


2 


61500 


Ni 


-0.12 ±0.11 


Q2348 


-01 c 


2 


61500 


Zn 


< 0.14 


Q2348 


- 1444 d 


2 


27900 


Cr 


+0.29 ±0.12 


Q2348 


- 1444 d 


2 


27900 


Ni 


< -0.28 


Q2348 


- 1444 d 


2 


27900 


Zn 


s; +0.28 



-2.01 ± 0.09 






1 


-2.01 ± 0.09 






1 


-2.01 ±0.09 






1 


-2.02 ±0.10 






3 


-2.20 ±0.13 






2 


-2.20 ±0.13 






2 


-2.20 ±0.13 






1 


-2.21 ±0.09 






< 3 


-2.21 ± 0.09 






1 


-2.21 ±0.09 






1 



"For the cases where two errors are given, they refer to random and systematic errors respectively. 
1 Molaro et aflfeOOOl). 

' ProchaskaetalJfeOOlh. 

' Dessauges-Zavadskv et alj i2007t) . 



in the middle panel of Fig.[4]as the open circles. In all cases, stel- 
lar abundances were derived from 1-D LTE models. Where neces- 
sary, the stellar measurements have been corrected to the abundance 
scale adopted in this work (see Section [3~Tl and TableO. 

Overlaid on the five panels of Fig. [4] are the DLA measure- 
ments and upper limits from Table [4] The red triangles are for 
DLAs with detections of Fe-peak element lines, or where a use- 
ful upper limit could be derived. The single blue triangle in each 
panel shows the mean value of [X/Fe] for the population of VMP 
DLAs without individual detections. The error bars on the values of 
[X/Fe] include both random and systematic uncertainties, and the 
ranges of [Fe/H] to which the stacks of undetected lines refer were 
estimated as explained in Section [3~4l We now comment separately 
on the individual Fe-peak element ratios. 



4.1 Titanium 

There are two resonance lines of Ti n accessible in the rest-frame ul- 
traviolet, at /11910.6123 and /U910.9538. The corresponding veloc- 
ity separation, Av = 53.7 km s is sufficient to avoid contamination 
between the two spectral regions over which we integrate, of width 
Av = ±25 km s centred on each line. In any case, any contamina- 
tion would have the result of making the upper limit we deduce for 



the Ti/Fe abundance more conservative. As can be seen from Fig. [3] 
(second row, left panel) even by stacking 16 Tin lines we do not 
achieve a detection; the 3<x upper limit we deduce, [Ti/Fe] < +0.16 
at [Fe/H] a* -2.2 ± 0.1 (Table® is plotted as a blue triangle in the 
top panel of Fig. [4] 

The derived limit on [Ti/Fe] is close to the lower envelope of 
[Ti/Fe] values in stars of comparable metallicity, derived from the 
analysis of Ti n linefl This does not leave much room for agree- 
ment between stars and DLAs, and the lack of any discernible ab- 
sorption beyond the lcr region near km s" 1 in the Ti/Fe panel 
of Fig. [3] further suggests some level of disagreement. Departures 
from LTE may result in sma ll corrections to th e stellar Ti/Fe ra- 
tios measured from Ti n lines jBergemannl201 lh . This may explain 
some of the mismatch between stars and DLAs. An alternative, 
or additional, possibility is that there may be a small residual de- 
gree of Ti depletion onto dust grains in VMP DLAs. Although it 
is usually assumed that DLAs with metallicities < 1/100 of so- 
lar ex hibit negligible depletion jPettini et alJI 19971 ; lAkerman et al .1 
l2005h . Ti is one of the most refractory elements, and is readily in- 
corporated into dust grains in the local interstellar medium (see e.g. 

2 We i gnore stellar abund ances derived from Ti i, as per the recommenda- 
tion by Bergemann ( 20 11). 
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Figure 4. Comparing Fe-peak element rat ios in VMP DLAs (triangles) and 
stars (open circles, s quares , and diamonds; lGratton et alj20O3llCavrel et alj 
[20041 and Lai et al.|| 2008| respectively). The open circles for [Co/Fe] are 
data from iBergemann. Pic kering. & Gehrer] j20ld) . The red triangles are 
for individual DLAs where the relevant absorption lines were detected, or 
where an informative upper limit could be obtained. The single blue triangle 
in each panel corresponds to the stack of DLAs in which no absorption by 
the element of interest was detected in individual systems. 



ISpitzer & Jenkinslll975h . We thus caution that Ti may suffer from a 
mild depletion in DLAs when [Fe/H] - -2.0. 



4.2 Chromium 

Chromium is accessible via the well known triplet 
Cr n ,L12056, 2062, 2066 (e.g. IPettini et al.lll99Lt) . In the present 
sample, we have four measures of [Cr/Fe] and two upper limits in 
individual DLAs, and a 5cr detection in the stack of DLAs without 



individual detections (see Table [4). These seven determinations 
are all in good mutual agreement, and they broadly agree with the 
stellar data, as can be seen in the second panel from the top in 
Fig. E] In the DLAs, [Cr/Fe] is approximately solar, or possibly 
slightly supersolar by ~ 0.1 dex, down to [Fe/H] ^ -2.5. 

Unfortunately, our survey does not extend to lower metal- 
licities where the stellar data appear to deviate from solar Cr/Fe, 
gradually decreasing to [Cr/Fe] =s -0.5 when [Fe/H] = -3.5 (when 
[Cr/Fe] is measured from Cri lines, as is the case for the data plot- 
ted in Fig. |4j- This decline may however be spurious, reflecting 
an increasing correction with decreasing metallicity for non-LTE 
(NLTE) effects to the Cri lines used in the analysis. On the basis 
of their NLTE calculations. iBergemann & Cescuttil d2010h have ar- 
gued that the [Cr/Fe] ratio remains close to solar down to the lowest 
metallicities. In future, it may be possible to verify this claim with 
observations of VMP DLAs of lower metallicity than those consid- 
ered here. This is an important check to perform because models 
of Population III nucleosynthesis interpret sub-solar Cr/Fe ratios 
as resulting from a deep er mass cut (i.e. larger explosion energy; 
lUmeda & Nomotoll2002l) . Such an explanation is in line with the 
observed rise in both the Co/Fe and Zn/Fe ratios at metallicities 
[Fe/H] < -2.5 (see below). 



4.3 Cobalt 

Together with Ti, Co is the least common of the Fe-peak elements 
considered here, with an abundance of only ~ 1 /350 of that of Fe 
in the Sun (Table [3). Thus, it is not surprising that we do not detect 
Co ii absorption in any of the 25 VMP DLAs in the present sample. 
This difficulty is offset to some degree by the availability of many 
Con transitions in the wavelength interval 1400-2000 A which is 
well observed in DLAs at z = 2-3 (see Table By stacking a 
total of 85 such features we are able to set an instructive 3<r upper 
limit of [Co/Fe] < +0.53 for the VMP DLA population at a mean 
metallicity [Fe/H] -2.2 (Table|4j. 

The earliest compilations of element ratios in the most metal- 
deficient stars in dicated that Co is enhanced re l ative to Fe when 
[Fe/H ] < -2.0 ( iRvan. Norris & Besselll 1 199ll: iMcWilliam et al.1 
[1991 . if the Co abundance is measured from Co i lines assuming 
LTE, as is the case for the data in Fig. [4] Even though later stud- 
ies have co nfirmed the trend o f increasing [Co /Fe] with decreasing 
metallicity dCavrel et al.l2004l:lLai et al.l2008l) . the true level of the 
enhancement of Co relative to Fe is still a matter of debate. The 
Co i lines appear to suffer from large (positive) corrections due to 
the neglect of modelling their hyperfine structure (HFS) as we ll 
as departures from LTE l lBergemann. Pickering. & Gehrenll2010h . 
In some cases, the offset is calculated to be as large as +0.6 dex. 
Applying these corrections would result in a substantial enhance- 
ment of Co, to a level that may be inconsistent with the upper limit 
to [Co/Fe] we have determined in the DLAs (see middle panel of 
Fig- EJ- Furthermore, model calculations of Co production in CC- 
SNe have difficulty in reproducing such extreme enhancements of 
the Co/Fe ratio. Some authors have interpreted the rise in [Co/Fe] 
at low metallicities as evi dence for high energy CCSN enrichment 
dUmeda & Nomotoll2002h . but even then there is still a shortfall in 
the model abundances compared to the stellar values corrected for 
NLTE and HFS effects. DLAs can provide new insights into this 
problem, and will undoubtedly play an important role in future to 
help settle the debate. 
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4.4 Nickel 

Ni is the second most abundant Fe-peak element in the Universe. 
This fact, in addition to the multitude of relatively strong transitions 
that are easily accessible (i.e. the lines are usually free from both 
the Lycf forest and telluric absorption), has resulted in Ni being the 
most commonly observed Fe-peak element in high redshift DLAs 
(besides Fe). For the sample of metal-poor DLAs considered here, 
we have seven detections and two upper limits in individual sys- 
tems, and a total of 52 lines contributing to the stack of undetected 
lines. 

In Galactic stars, the relative abundances of Ni and Fe remain 
essentially constant over the entire range of metallicity probed. Al- 
though there are notable exceptions to the rule, the typical stellar 
[Ni/Fe] value is solar ([(Ni/Fe>]= -0.01), in very good agreement 
with the value [Ni/Fe] = -0.03 ± 0.04 ± 0.03 we derive for VMP 
DLAs at a mean [Fe/H] ^ -2.2 (see Fig. [4). It is interesting, how- 
ever, that most of the systems with detected Ni n absorption exhibit 
slightly sub-solar ratios, except for one DLA with an unusually 
high Ni/Fe ratio. It appears that both stars and DLAs show signif- 
icant variation in their Ni/Fe ratios, and DLAs preferentially have 
sub-solar values of [Ni/Fe]. 



4.5 Zinc 

Zn is rarely detected in DLAs when the metallicity is less than 
1/100 of solar; w e have only one de tected system in our sam- 
ple (reported by iMolaro et alj |2000i) in a DLA with one of 
the largest recorded values of neutral hydrogen column density, 
log N(H i)/cnT 2 = 21.41. In addition, there are a further two sys- 
tems in the literature with useful upper limits. From the stacked 
spectrum of 1 1 unblended Zn n lines in DLAs without individual 
detections, we deduce a 3cr upper limit [Zn/Fe] < -0.04 at a mean 
[Fe/H] -2.3. 

While [Zn/Fe] is clo se to solar in stars with [Fe/H] > -2.0 
dNissen et alj|2004 120071) . it rises to supersolar values at lower 
m etallicities (see botto m panel of Fig. |4j. T his trend, firs t noted 
bv lPrimas etail ( feOOOh and later confirmed bv lCavrel et~aH J2004h 
with a larger sample of halo stars, is unlikely to be due to the de- 
tails of the stellar atmosph ere modelling. The curren tly favoured 
interpretation, suggested bv lUmeda & Nomotd d2002l) , is that it re- 
flects nucleosynthesis by hypernovae — supemovae with an order 
of magnitude larger explosion energy than normal. The higher en- 
ergy explosion results in a stronger o--rich freezeout, which in turn 
favours the production of Zn. 

Although the statistics for VMP DLAs are limited, our 3cr up- 
per limit [Zn/Fe] < -0.04 for the population of DLAs with a mean 
[Fe/H] =2 -2.3 (blue triangle in the bottom panel of Fig. [4]( is not 
consistent with an enhancement of Zn at these metallicities. There 
is no hint of Zn n absorption in the stack shown in the bottom panel 
of Fig. [3] suggesting that the true value of [Zn/Fe] in VMP DLAs 
may well be sub-solar and thus different from the values found in 
stars of comparable metallicity. Given that Zn has the potential to 
be a powerful probe of the explosion energy, it is now important 
to increase the sample of DLAs with measurable [Zn/Fe] down to 
at least [Fe/H] = -3.0 for a more instructive comparison with the 
stellar data than is possible at the moment. 



5 EXPLOSION ENERGY OF TYPE-II SUPERNOVAE 

It was shown in ICookeetal.N2011bh that the chemical composi- 
tion of the most metal-poor DLAs is consistent with the hypothesis 
that they were enriched by massive stars that ended their lives as 
core-collapse supernovae. In this section, we will attempt to esti- 
mate the typical explosion energy of these supernovae by compar- 
ing the abundances of Fe-peak elements deduced above with the 
relative yields calculated with models of explosive nucleosynthesis 
in metal-free stars. It has been known for over a decade that the 
relative abundance of the Fe-peak elements at low metallicity can 
provide clues to the physic s of a core-collapse supern ova explosion 
dWooslev & Weave il l 1995k lUmeda & Nomotd |2002|) . Progress in 
this area is highly desirable, because we still lack a physical under- 
standing of the me chanism that drives the explosion, despite many 
recen t efforts (e.g. lMe zzacappa 2005; Ott. O'Connor. & Dasguptal 
|20T 1). In the current generation of massive star nucleosynthesis 
models this shortcoming is circumvented by parameterising the 
physics of the explosion in order to calculate the resulting element 
yields. The parameterisation is typically achieved by employing a 
simple description of the mixing that occurs between the stellar lay- 
ers (due to either Rayleigh-Taylor or rotationally-induced mixing), 
and specifying a mass-coordinate (or mass-cut) that defines the ma- 
terial that either escapes the binding energy of the exploding star, 
or falls back onto the newly formed compact object. 

Perhaps the most d etailed CCSN exp l osive nucleosynthesis 
calculations are those by iHeger & Wooslevl d2010h . These authors 
have investigated the explosive nucleosynthesis in massive stars 
(M = 10-100 M Q ) with a fine mass-resolution (> 0.1 M G ); the 
entire suite consists of 120 stellar models. Each star undergoes a 
supernova explosion, which is simulated as a moving piston that de- 
posits momentum at a specified mass-coordinate. We have adopted 
the standard case which sets the piston to be at the base of the 
oxygen burning shell, corresponding to an entropy per baryon of 
- 4k B . Each exploding star is simulated with 10 different values of 
the final kinetic energy at infinity, taken to be in the range (0.3 - 10) 
xlO 51 erg. A simple prescription of the mixing that occurs between 
the stellar layers during the explosion is also implemented. This is 
achieved by running a boxcar filter with a variety of widths (14 in 
total), ranging from no mixing at all to the maximum filter width 
of 25% of the He core size. The entire suite of models consists of a 
total of 16 800 combinations of these parameters. 

By enforcing the selection criterion that the ejected mass of 
56 Ni is wit hin the observed range of 0.001 < M( 56 Ni) < 1 M G (see 
figure 1 of iTominaga. Umeda. & Nomotdl2007l) . we are left with 
a subset of 10790 'realistic' models. We subdivide these models 
into four arbitrary groups depending on the final kinetic energy 
at infinity: (1) the hypernovae, with energies in the range (5.0- 
10.0) x 10 51 erg; (2) the high-energy supernovae, with energies in 
the range (1.8-3.0) x 10 51 erg; (3) the typical supernovae, with en- 
ergies in the range (0.9-1. 5)xl0 51 erg; and (4) the faint supernovae, 
with energies in the range (0.3-0.6) x 10 51 erg. In Fig. [5] these four 
groups are colour-coded in red, green, blue, and black respectively. 
Given the large number of models considered, in Fig. [5] we repre- 
sent each of these four explosion energy ranges with two ellipses 
that enclose, respectively, 68% (solid lines) and 95% (dotted lines) 
of the models in that range. In this way, we can learn how likely a 
certain explosion energy is, given two measured Fe-peak element 
ratios. For all panels in the figure, we also show the ranges of Fe- 
peak element ratios determined from the DLA stacks (displayed as 
shaded boxes). 

As can be seen from the top panels in Fig. [5] the abundances 
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Figure 5. Comparison of the explosive nucleosynthesis calculations bv lHeger&Wooslevl l2O10t> with elements ratios measured in VMP DLAs. In each 
two-ratio plane, the theoretical values are shown with ellipses enclosing 68% (solid ellipses) and 95% (dotted ellipses) of the models considered in each band 
of supernova explosion energy. Four bands are considered, colour-coded red, green, blue, and black. Broadly speaking, red corresponds to hypernovae, green 
is for high energy Type II supernovae, blue represents a typical Type II supernova, and black corresponds to the faint supernovae (see text for further details 
and the specific cuts we have employed). The blue shaded region in each panel shows the values measured in the stack of VMP DLAs. 



of Ti, Cr, and Ni relative to iron all have very little discrimina- 
tory power on the explosion energy (the ellipses largely overlap 
on the planes defined by the ratios of these elements). However, 
it is at least reassuring that model calculations and VMP DLA 
measurements of the Cr/Fe and Ni/Fe ratios are in good agree- 
ment with each other. The Co/Fe and Zn/Fe ratios are somewhat 
more sensitive to the explosion energy; the latter appears to be the 
better discriminant, and to illustrate this point better we show in 
Fig. [6] the full di s tributi on of [Zn/Fe] values expected from the 
iHeaer & Wooslevl l l2010h models with the same colour-coding as 
in Fig. [5] It is also interesti ng to note in Fig.[5] t hat Ni appears to be 
co-produced with Co in the lHeger & Wooslevl j201Ct) models. 

The relationship bet ween Zn/Fe and ex p losion energy has 
been previously noted by lUmeda & Nomotol j2002l) . These au- 
thors invoked the increased importance of hypernovae with de- 
creasing metallicity to explain the gradual rise to [Zn/Fe] = +0.7 



at [Fe/H] ~ -4.0 discovered in the m ost metal-poor halo stars 
jPrimas et alj|2000l : ICavrel et al.ll2004h. However, ev en the most 
energetic supernovae in the iHeger & Wooslevl < l2oToh models un- 
derproduce ZjT_rompared_to_the enhan ced values reported by 
ICavrel et all d2004h andlLai et al.N2008t) . To account for this de- 
ficiency, Heger & Wooslevl J2010h caution that Zn can also be 
produced by other (unaccounted for) sources, such as the neu- 
trino wind from a newly formed neutron star (i.e. the 'hot bub - 
ble' component; f or a r elevant discussion, see iPruet et al. IIH1). 
Ilzutani & Umedd d2010h . however, point out that such a scenario 
requires a contrived degree of fine-tuning for the explosion param- 
eters and favour hypernovae as the more natural explanation for the 
elevated Zn/Fe ratios. 

Thus, the jury is still out regarding the origin of the high Zn/Fe 
ratios measured in the most-metal p oor stars. As has be en shown 
to be the case for other elements (e.g. lCooke et al .1201 lbl) . comple- 
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Figure 6. The distri bution of [Zn/Fe] values in the subset of 
iHeger & Wooslevl feOld) models considered. The four sets of models are 
colour-coded according to the kinetic energy released by the supemovae. 



mentary observations in VMP DLAs have the potential of resolving 
the impasse. Unfortunately, the Zn n lines in VMP DLAs are very 
weak, and their detection with current instrumentation will require 
the rare combination of bright background QSOs and unusually 
high column densities of neutral hydrogen. We estimate that DLAs 
with A'(Hi) > 10 21 cnr 2 are required to have a chance of detecting 
the Znn lines in the metallicity regime [Fe/H] < -2.5; such DLAs 
do exist, but all known cases are located in front of QSOs too faint 
for their spectra to be recorded at the required S/N and resolution 
with current instrumentation. However, measuring the Zn/Fe ratio 
in the most metal-poor DLAs is a goal that can be readily achieved 
with the forthcoming generation of 30-m class telescopes. This a 
very exciting prospect. Only by measuring the relative abundance 
of the Fe-peak elements in DLAs with [Fe/H] < -3.0, can we con- 
firm the trend of [Zn/Fe] at the lowest metallicities, and thereby pin 
down the typical energy released by the CCSNe of the first stars. 

In the more immediate future, it may be possible to identify 
gas clouds that were enriched by hypernovae, or high-energy su- 
pernovae, by searching for D LAs with a modest [Ni/Fe] enhance- 
ment. The only models in the IHeger & Wooslevl feOld) suite that 
are able to reproduce even modest super-solar values of [Ni/Fe] 
are the hypernovae. A system with such properties may be present 
in our sample: the DLA at z abs = 2.84430 towards J1944+7705 
has [Ni/Fe] = +0.25 and is clearly an outlier relative to the typi- 
cal DLA population in this metallicity regime (see Fig.Q- Deeper 
data are required to test if this system also exhibits enhanced Zn/Fe 
and Co/Fe ratios, which are the only sure diagnostics that can 
confirm such a scenario; at present we place 3cr upper limits of 
[Zn/Fe] < +0.39 and [Co/Fe] < +1.06. 

The primary motivation of this work was to estimate the en- 
ergy released by the stars that enriched the most metal-poor DLAs. 
The results presented here are only a first step in this direction and 
have mostly highlighted the need for further observations of DLAs 
at even lower metallicities than the regime probed so far. Based on 
the low value of [Zn/Fe] we have found, we tentatively conclude 
that VMP DLAs were probably enriched by massive stars that re- 
leased < 3.0 x 10 51 erg of energy, However, we caution that this 
conclusion is based on the assumption that these absorbers were en- 
riched only by Population III stars (and that the nucleosynthesis by 
such stars is accurately reproduced by current models). It remains 
to be tested how the relative Zn and Fe yields depend on the metal- 
licity of the progenitor star, as well as the explosion energy. This is 



a goal for future, high mass-resolution, metallicity-dependent stel- 
lar yield calculations. 



6 SUMMARY AND CONCLUSIONS 

We have presented the first dedicated survey of Fe-peak elements in 
very metal-poor DLAs at redshifts z ~ 2-3, focussing in particular 
on the ratios of Ti, Cr, Co, Ni, Zn to Fe. Based on the analysis of 
the relative abundances of these elements in a sample of 25 DLAs 
with [Fe/H] < 2.0, we draw the following conclusions: 

(i) In nine of the 25 DLAs, we could measure element ratios, or 
useful upper limits, directly from detected absorption lines. In the 
majority of DLAs, however, the lines are too weak to be detected in 
available spectra. To deal with such cases, we have devised a form 
of spectral stacking that allows us to measure the Fe-peak element 
ratios typical of the VMP DLA population as a whole, at mean 
metallicities [Fe/H] =s -2.2 to -2.3. 

(ii) Comparing the Fe-peak element ratios deduced in VMP DLAs 
with the values measured in metal-poor stars of the Galactic halo, 
we found fair agreement for [Cr/Fe] and [Ni/Fe]. Furthermore, the 
upper limit on [Co/Fe] in VMP DLAs is consistent with the stellar 
values, provided the non-LTE corrections to the latter are not too 
severe. On the other hand, the upper limits we place on the [Ti/Fe] 
and [Zn/Fe] values in VMP DLAs are barely consistent with the 
stellar measurements; the DLA upper limits lie near the lower en- 
velope of these ratios in halo stars. Whether there is indeed a dis- 
crepancy or not can only be confirmed by future observations of 
DLAs with [Fe/H] < -2.5 , this being the metallicity regime where 
the stellar trends exhibit marked deviations from the solar ratios. 

(iii) By comparing the Fe-peak element ratios in the DLAs 
to the values computed with the most detailed current models 
of nucleosynthesis by metal-free stars, we tentatively concluded 
that the DLAs in our sample were probably enriched by massive 
stars that released < 3.0 x 10 51 erg when they exploded as core- 
collapse supernovae. An extension of models such as those by 
lHeger& Wooslevl fcO id) to Population II stars is required to test if 
the Fe-peak element yields, particularly those of Zn and Fe, depend 
on metallicity as well as the stellar mass and explosion parameters. 

(iv) We have uncovered a relatively Ni-rich DLA, with [Ni/Fe] = 
+0.25, a result which may hint at an origin of these elements in 
hypernovae. Such an interpretation is consistent with the 3<r upper 
limits [Co/Fe] < +1.06 and [Zn/Fe] < +0.39 we can set for this 
system with available data. 

The present belief is that the most metal-poor Galactic stars in- 
corporated the metals synthesised by an earlier generation of mas- 
sive stars that had ended their lives as hypernovae; however, we 
have found little evidence so far for this hypothesis in the chemi- 
cal composition of the most metal-poor DLAs. Future observations 
with the forthcoming generation of 30-m class telescopes, equipped 
with echelle spectrographs, will allow us to confidently measure the 
relative abundances of Fe-peak elements in individual DLAs with 
[Fe/H] < 1 / 1000 of solar. We suspect that only then will we be able 
to reliably determine the energy released by the first stars. 
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